Core phase structure of cosmic strings and monopoles 



M. Axenides 

Institute of Nuclear Physics, N.C.S.R. Demokritos, 153 10, Athens; Greece 
e-mail: axenides@gr3801.nrcps.ariadne-t.gr 

L. Perivolaropoulos and T.N. Tomaras 
Department of Physics and Institute of Plasma Physics, University of Crete and FO.R.T.H. 
P.O.Box 2208, 710 03 Heraklion, Crete; Greece 
e-mail: leandros@physics.uch.gr, tomaras@physics.uch.gr 
(February 1, 2008) 



00 
OS 



in 



> 
in 
in 
m 
m : 
o . 
oo ■ 

5 : 

a: 
i 

^ ■ 



X 



Global and local symmetries may or may not be restored 
inside topological defects depending upon the values of the 
parameters of the model. A detailed study of this parameter 
dependence of the core structure of strings and monopoles is 
presented in the context of simple models. 



I. INTRODUCTION 



A lot of theoretical and experimental effort has been 
devoted over the past twenty years to the search for de- 
fects predicted in field theory models and to the study 
of their dynamics and role in particle physics and cos- 
mology pj] |Q. Despite of the fact that todate there is 
no direct evidence of any such object, their potential im- 
portance in cosmology (structure formation, baryogene- 
sis and the interest in the non-perturbative structure 
of field theories has provided enough motivation for the- 
orists to search for all possible extended solutions and 
to analyze their physical properties. The absolutely sta- 
ble topological solitons were studied first ||, and were 
followed by the discovery of the metastable semilocal de- 
fects p|, j^], the electroweak strings ||, and the ribbons 
in realistic models of high energy physics, including 
topologically trivial ones. 

The role of solitons is known to depend to some ex- 
tent on the detailed profile and the unbroken symmetries 
in their interior. For instance, the cosmological evolution 
of a superconducting string Q differs significantly from 
that of a normal string . The purpose of this note is to 
study the parameter dependence of the core structure of 
strings and monopoles. As a first step, we restrict our- 
selves to models which are simple enough to be analyzed 
in detail, and which furthermore could easily be embed- 
ded into larger realistic theories. This analysis is par- 
ticularly interesting in the context of Cosmology, where 
due to the temperature dependence of the parameters of 
field theories, one may encounter during the cosmologi- 
cal evolution " phase transitions" inside the cores of topo- 
logical defects. In fact, laboratory experiments on 3 He, 
designed to investigate the physics of phase transitions 



in the Early Universe, have explicitly provided us with 
strong experimental evidence for defect-core transitions 
in the interior of vortices which appear in the superfluid 
3 He — B phase fuj. 

In section II we study the core structure of cosmic 
strings. The analysis is done in the context of a U(l) 
gauge model with two classically relevant parameters and 
with the original semilocal string as a limiting case. A 
similar analysis is performed in section III where the 't 
Hooft-Polyakov monopole is embedded in an SU(2) gauge 
model with an extended Higgs sector. The parameter 
space is divided in regions corresponding to the two pos- 
sibilities of the monopole-core global symmetry. Possible 
extensions and applications are commented upon in the 
final discussion section. 



II. EMBEDDED NIELSEN-OLESEN VORTEX 

Consider the simple extension of the Abelian Higgs 
model with two equally charged scalars $1 and $2 de- 
scribed by the lagrangian density 
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( J D AI $) t ^$-U($i,$ 2 



(1) 



where $ = ($i,$ 2 ), B 
igB^, and with T3 the third Pauli matrix 



= d^B v — d u B^, 



D, 



= du- 



M 2 to 2 A 

U($i,$ 2 ) = ^-$ t $-— $ t r 3 $ + -($ t $) 2 (2) 



For generic values of the parameters the vacuum has 
= (Af 2 + to 2 )/ A and $ 2 = 0. The model is an 
example of a U(l) gauge theory with the gauge group 
broken spontaneously to the identity, and as such it sup- 
ports the existence of topologically stable strings. $1, 
the field with non-vanishing vacuum expectation value, 
carries their winding and vanishes at the center. As for 
7 it vanishes at infinity, but its profile inside the string 
has no geometric or topological constraints and is deter- 
mined dynamically by the field equations. In this section 
we will study in detail the structure of the string core 
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and how it varies with the parameters of the model. One 
might imagine (|l|) , (||) embedded in a more realistic the- 
ory with $2 coupled to electromagnetism. A non-zero 
core value of $2 in this case would render the string su- 
perconducting with well known observable effects S. 

Of the four parameters of the model, one sets the 
scale, another may be pulled outside of the action to play 
the role of the semiclassical parameter Ti and there remain 
two classically relevant ones. Specifically, by rescaling 
fields and coordinates according to 

$ 1(2) -> M$ 1(2) /VA, -> MBjg, Xl -> xt/M (3) 



one is left with the parameters 



^- and (3 = —=. 



(4) 



Various limiting cases of ([!]), (|^) have been studied be- 
fore. First, for both a and j3 equal to zero a simple 
scaling argument shows that the model does not possess 
any kind of stable defect solution. Second, for a = one 
obtains an SU(2) 9 ; (, a ; x U(l)z oca ; symmetric model, with 
the U(l) gauge field distinguishing an S 1 fiber out of the 
S 3 vacuum. The "strength" of this fibration is propor- 
tional to (3. For (3 > (3 = \[2 the embedded Nielsen- 
Olesen configuration i.e. with $2 = is a classically 
stable solution of the model jll|] , and is the first example 
of a semilocal vortex studied before [Q . For smaller val- 
ues of (3 an instability arises due to the development of 
a non-vanishing $2 inside the core and the string blows 
up to the vacuum. Finally, the global model with (3 = 
was studied recently in ]l2|] . In analogy to the U(l) gaug- 
ing of the previous limiting case, one may think of the 
term |<E>i| 2 — \$2\ 2 in (||) as defining an S 1 "scalar" fibra- 
tion with "strength" a on the S 3 vacuum manifold of the 
a = model. Topologically stable global strings (with 
logarithmically infinite energy per unit length) exist in 
this case for all positive values of a. Furthermore, for 
a > cto ~ 0.4 they have $2 = 0, while $2^0 when 
a < ctQ. 

The simplest guess for the parameter dependence of 
the core structure of the strings in (|l|), (||) is then the fol- 
lowing: In the (a, (3) plane there is a curve connecting the 
critical points {a®, 0) and (0, (3$) on the two axes, outside 
of which the embedded Nielsen-Olesen strings are stable, 
and inside of which the stable strings are characterized 
by non- vanishing $2. This picture will indeed be con- 
firmed below and the curve a cr it((3) will be determined 
numerically. 

The ansatz for the string solutions with winding num- 
ber n has the axially symmetric form 



<f 1 =f(p)e in0 , $2 = G(p), B 



B(p) 



(5) 



and the corresponding energy functional is in units of 
M i /2X 



J + u +^2~3" H -3. J 



f3 2 p 2 



E = J d 2 x 

+ K g2 + \(f + g 2 ) 2 - (i + a2 )f - (1 - a2 ) G2 



(6) 



It leads to the following field equations: 



f" 



-f+(l + a 2 -f 2 -G 2 )f = 



B' 



B" - — + f3 2 f(n~ B)~ (3 2 G 2 B = Q (7) 
P 

G" + — -^G+(l-a 2 -f 2 - G 2 )G = 
P P 

We restrict our analysis to the minimal n — 1 case. In the 
end we comment briefly on the results for higher values 
of n. Finiteness of the energy and the field equations at 
the origin imply the boundary conditions 



/(oo) = yja 2 + 1, G(oo) = 0, B(oo) 
/(0) = 0, G'(0) = 0, B(0) = 



(8) 



For all values of the parameters a and f3, equations 
(0) admit the embedded Nielsen-Olesen solution with 
G(p) = 0. As explained above, this is known to be clas- 
sically stable on the a-axis for a > ao and on the /3-axis 
for (3 > (3q. For arbitrary (a, (3), the region of stability 
of the corresponding solution is determined by the re- 
quirement that in the expansion of the energy functional 
around it the first non-trivial term in SE is strictly pos- 
itive for all field variations. To quadratic order in field 
changes SE = SE^o + SEq, where SEmo is independent 
of G(p) and identical in form with the perturbation ob- 
tained for the topologically stable Nielsen-Olesen vortex. 
Thus SEno > for all Sf and 8 A and may be ignored. 
SEg on the other hand represents the quadratic correc- 
tion to the energy of the Nielsen-Olesen string due to the 
excitation of G(p). It may be written in the form: 



SE G = 



d 2 x GOG 



with 



O 



1 d d B 2 „ 2 2 
-— p— + — + f 2 + a 2 
p dp dp p A 



(9) 



(10) 



where /, B are the Nielsen-Olesen fields obtained by solv- 
ing the first two equations of system (]?]) with G — 0. For 
stability of the embedded Nielsen-Olesen vortex we re- 
quire that SEg > for arbitrary small perturbation G. 
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This is equivalent to demanding that the operator O has 
no negative eigenvalues. The region of stability of the 
embedded Nielsen-Olesen vortex was determined numer- 
ically. The method is straightforward. We pick a point 
in the (a, (3) plane, set G = and solve the remaining 
equations (^) to find the corresponding Nielsen-Olesen 
solution. We then construct the operator O and consider 
the eigenvalue problem 



OG = lu 2 G 



(11) 
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Figure 1: The string-core structure for all values of 
parameters a, (3. Inside the region bounded by the solid 
line $2 7^ 0, while outside it $2 vanishes. 

We use a fourth order Runge-Kutta routine with ini- 
tial conditions G(0) = 1, G'(0) = to investigate if there 
is a bound state in equation (|Tl]). If there is a po such 
that for p > po we have G(p) < with the above initial 
condition, then there is clearly a bound state for the con- 
sidered values of a, (3. This implies an instability of the 
embedded Nielsen-Olesen vortex towards another stable 
configuration with non-zero order parameter in the core. 
In Figure 1 we plot the parameter space (a, (3) and dis- 
play the regions of stability and instability of the embed- 
ded Nielsen-Olesen vortex. We have repeated the above 
analysis to the \n\ = 2 case. As shown in Figure 1 in- 
creasing \n\ leads to an expansion of the superconducting 
region. This is related to the fact that the behaviour of 
/ near zero is / ~ r> n >. Thus increasing \n\ leads to a 
broadening of the potential of the Schroedinger operator 
O and favours the existence of negative eigenvalues. 

Figures 2 and 3 contain the results of the numerical 
integration via a relaxation method of the system 



In Figure 2 we show a stable embedded Nielsen-Olesen 
configuration and in Figure 3 an example of a vortex solu- 
tion for parameter values in the superconducting region. 
The latter has a non-zero order parameter in the core. 
For ci / this relaxed configuration is stable due to a 
non-trivial topology. For a — there is no topological 
stability and the Nielsen-Olesen configuration spreads its 
flux to infinity pT|. 
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Figure 2: For a = 0.1 and (3 — 1.3 the stable string 
configuration shown here is the embedded Nielsen-Olesen 
solution. 



1.2 



1.0 
0.8 
0.7 
0.5 
0.3 
0.2 
0.0 







f(p) 




V '' 


B(p) 




a=0.1 




(3=0.6 




("'( P) 



0.0 5.0 10.0 15.0 20.0 25.0 30.0 
P 



Figure 3: The profile of the stable string for a = 0.1 
and (3 = 0.6. It has $2 7^ inside the core. 
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To leading order in the temperature T the parameter 
a becomes T-dependent a 2 (T) = -m 2 /[-M 2 + (2A + 
g 2 )T 2 /4], while 8 is constant. Strings form at T slightly 
below the critical temperature Tq. a 2 {T ~ Tc) is very 
large and the strings start off normal. When T is such 
that a(T) = a cr it(8) for the given 8 the string becomes 
superconducting. 



III. STABILITY OF THE EMBEDDED 
MONOPOLE 

It is straightforward to apply the above method to 
the study of the fate of global symmetries inside the 
monopole core. We choose to work in the context of 
the simple model 

C = -~F° u F a,i " + ^D^ a D^^ a + i(a M $ 4 ) 2 - V (12) 

describing the dynamics of an 0(3) gauge field A a ^ cou- 
pled to the scalar triplet <& a a=l,2,3, which in addition 
interacts with the gauge singlet <I> 4 through the potential 

V = ^($ a $ a + $ 4 $ 4 - v 2 ) 2 - —($ a $ a - $ 4 $ 4 ) (13) 



the string case can only be decided by solving the field 
equations, the symmetry inside the monopole core will 
be either 0(3) (G(r) ^ 0) or the full 0(4) (G(r) = 0). 
Which one of the two is realized depends on the values of 
the parameters of the model. Using the same numerical 
method as in the string case we give below a complete 
map of the parameter space based on the core symmetry 
of the corresponding magnetic monopole. 

It is convenient to define K(r) = l — 8rW(r), in which 
case the field equations for the three unknown functions 
of the ansatz take the form: 

/" + v " % r2 + (1 + &2 - /2 ~ ° 2)f = 

K" - K[K l~ l) - ~8 2 fK = (16) 
or" 

G" + — + (1 - a 2 - f - G 2 )G = 
r 

while the corresponding boundary conditions, dictated 
by the finiteness of the energy and the field equations at 
the origin, are 

/(oo) = v/l + « 2 , G(oo) = 0, ^(co)=0 
/(0)=0, G'(0)=0, K{0) = 1 



The field strength and the covariant derivative are given 



by F- 



le abc A b A c and D ff,a = Q S>a. 



ge abc A b ^ c , respectively. 

The above is a simple extension of the Georgi- 
Glashow 0(3) model with the two classically relevant 
parameters 



a 



m 
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and 8 = —= 



(14) 



as revealed after the rescalling x l — » x l /(v^/X), <i> a 



V<P 4 



and Af, 



For a = (3 — ((lj) possesses an 0(4) global symme- 
try. A non-zero m breaks this symmetry explicitly to the 
gauged 0(3) subgroup, and this in turn is spontaneously 
broken to 0(2) by the vacuum of the model. According 
to the standard lore the model admits for generic values 
of a and 8 a whole tower of topologically stable magnetic 
monopoles. The minimal one in its "ground state" has 
the spherically symmetric form 



$ a = S la -f(r), Af=e aij -W(r 



(15) 



* 4 



G(r) 



with f(r) and W(r) necessarily vanishing at the origin 
r = 0. Depending on the profile of G(r), which as in 
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Figure 4: The dependence on the parameters a and 
8 of the monopole-core global symmetry. 

We now look for instability modes of the embedded 
't Hooft-Polyakov monopole solution (i.e. with G(r) = 
at all r) of equations (|l6|). Using the same approach as 
in the case of the Nielsen-Olesen vortex we obtain the 
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linearized eigenvalue problem corresponding to the last 
equation in ( |l6| ) 

- G" - —— + (/ 2 - 1 + d?)G = uj 2 G (18) 
r 

where / is obtained from the 't Hooft-Polyakov monopole 
solution i.e. solving the system of the first two equations 
in ( |l6| ) with G = 0. Notice that contrary to the string 
case there is no coupling of the gauge field to G. In mod- 
els where such coupling exists and the gauge symmetry 
is 0(4) the embedded 't Hooft-Polyakov monopole is al- 
ways unstable due to the Brandt-Neri-Coleman instabil- 
ity Jl3| . This instability is realized as a screening of the 
long range magnetic field of the monopole by the other 
massless gauge fields of the theory 14 1. No such insta- 
bility exists in the model under discussion because the 
gauge symmetry is 0(3) and not 0(4). 
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Figure 5: For a = 0.2 and j3 = 0.6 the stable 
monopole is the embedded 't Hooft-Polyakov solution 
with G(r)=0. 



We have solved the eigenvalue problem (|18|) for sev- 
eral parameter pairs (a, (3). The parameter space shown 
in Figure 4 is divided into two regions according to the 
symmetry properties of the core of the stable monopole 
solution. For strong scalar and gauge fibrations i.e. for 
a, /3 large, the core symmetry of the stable monopole so- 
lution is 0(4). Negative modes (instabilities) towards 
an 0(3) symmetric core develop for weak fibrations as 
in the case of the embedded Nielsen-Olesen vortex. The 
embedded 't Hooft-Polyakov monopole in the semilocal 
limit was first discussed in ||. Our results confirm the 
qualitative analysis presented there. 



We solved numerically equations (16) and (17) using 
a relaxation method for a variety of parameter values. 
In Figure 5 an example of a stable embedded 't Hooft- 
Polyakov monopole solution is shown, while in Figure 6 
we plot the profile of a stable monopole solution with 
0(3) symmetric core. 
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Figure 6: The stable monopole shown here for a = 0.2 
and f3 = 0.1 has $ 4 excited inside the core. 



IV. CONCLUSION 

The dependence of the core phase structure of flux 
vortices and magnetic monopoles was studied explicitly 
in the context of simple gauge models carrying these 
topological defects. The combination of "scalar" and 
"gauge" fibrations on the vacuum manifold leads to in- 
stabilities of the core of the embedded Nielsen-Olesen 
vortices and 't Hooft-Polyakov monopoles, which result 
in a non-trivial interior in the stable solution. The mod- 
els discussed here are very simple and allow for only two 
phases in the corresponding defect interiors. A larger 
variety of core phases is of course expected in realistic 
models with more Higgs multiplets and richer symmetry 
pattern. The cosmological implications of the transitions 
between these different core possibilities deserve further 
study. 
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